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UNIQUENESS OF THE COMPACTLY SUPPORTED WEAK 
SOLUTIONS OF THE RELATIVISTIC VLASOV-DARWIN 

SYSTEM 

REINEL SOSPEDRA-ALFONSO AND MARTIAL AGUEH 



Abstract. We use optimal transportation techniques to show uniqueness of 
the compactly supported weak solutions of the relativistic Vlasov-Darwin sys- 
tem. Our proof extends the method used by Loeper in [8] to obtain uniqueness 
results for the Vlasov-Poisson system. 



1. Introduction 



r^ I The relativistic Vlasov-Darwin (RVD) system describes the evolution of a col- 

lisionless plasma whose particles interact only through the electromagnetic field 
they induce. In contrast to the Vlasov-Maxwell system, the particle interaction is 
assumed to be a low-order relativistic correction (i.e., the Darwin approximation) 
of the full Maxwell case. In [T3], we showed that the model equations for a single- 
species can be written in terms of a scalar and vector potentials {^,A) according 
to 



> 

(N 

00 



(1.1) 5J + v^-V,/-[V$-c-iv^VA']-Vp/ = 0, VA P c A 



l + c-2 \p-c^^A\' 



coupled with 






(1.2) 


^{t,x) 


Jr^ \y - x\ 


(1.3) 


A{t,x) 


2cJm3 \y-x\ 


via 






(1.4) 


p{t,x)^ 1 


f{t,x,p)dp, JA{t,x)^\ VAf(t,x,p)dp 



Here / = /(i, a;,p) denotes the one-particle distribution function at a time t e]0, cxd[, 
position a; £ M'^ and (generalized) momentum p G M'^. The scalar and vector 
potentials $ and A are induced by / via the charge and current densities p and 
JA, respectively. The relativistic particle velocity is denoted by va, and satisfies 
\va\ < c where c is the speed of light. For simplicity, we have set the charge and 
mass of the particles to one. Also, we have denoted the 3-by-3 identity matrix as 
id and the unit vector uj = {y ~ x) / \y — x\. As usual, repeated indexes means 
summation, which is the case in (|l.ip . and the symbol O in (|1.3p stands for the 
tensor product. We define the weak solutions to the RVD system as follows: 
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2 REINEL SOSPEDRA-ALFONSO AND MARTIAL AGUEH 

Definition 1. Let fo e L^ n L°°{R^;R) and /o > 0. For T > consider f G 
C([0,T[,i°°(R6) -^ w*;M.) such that / > 0. We call f a weak solution of the 
relativistic Vlasov- Darwin system with initial datum /o, if 

I I fit,x,p)dxdp^ / / fo{x,p)dxdp Vie[0,r[, 

$ and A are given by mM)-(l^, and for all tp G C^{[0,T[xR^;R) we have 

I I I f{dtip + VA-V^(p-[\/^-c-^v\\/A']-\/p^}{t,x,p)dtdxdp 

Jo JR3 JR3 



(1.5) =-/ / fo{x,p)ip{0,x,p)dxdp. 

JR3 JR3 

The notation (p € C^([0,T[xR6;M) means that ip is of class C°°, ip{T) = and 
ip{t) has a compact support in R^ for all t G [0, T[. Also, L°°(M^) — -u;* means that 
the space L°°(R^) is equipped with the weak-* topology. 

The global in time existence of weak solutions to the RVD system without size 
restriction on the Cauchy datum was established by Pallard in [12] . This result uses 
the formulation of the RVD system in terms of the electromagnetic field, which is 
formally equivalent to (jl.ip - (|1.4p : see details in [14]. Here we do not discuss the 
existence problem nor the equivalence of the two formulations in the context of 
weak solutions. We are only concerned with the uniqueness of solutions to (jl.ip - 
(|1.4p in the sense of Definition [TJ We emphasize, however, that for a bounded 
charge density (and in particular for a compactly supported distribution function) , 
Definition [1] makes perfect sense. We postpone this discussion to the next section 
(see Lemmas [2] and HI . 

If we formally let c ^^ cxd in (|l.ip - (ll.4l) . the model equations reduce to the Vlasov- 
Poisson system, which is the zeroth-order relativistic correction of the Vlasov- 
Maxwell system. A uniqueness result for weak solutions of the Vlasov-Poisson 
system based on optimal transportation, under the assumption that the charge 
density remains bounded, was established by Loeper in [8], (see also Robert [13] 
for the uniqueness of compactly supported weak solutions). Here, we extend the 
method of [Sj to produce a uniqueness of weak solutions to the RVD system under 
the assumption of a compactly supported distribution function. To the best of our 
knowledge, this is the first uniqueness result for weak solutions of the RVD system. 
Our main result is as follows: 

Theorem 1. Let foGL^n L°°{R^;M.), /o > 0, have compact support and let f be 
a weak solution on [0,T[ of the RVD system il.l]) - (L^ in the sense of Definition 
[7J such that /It^Q = /o- If for all t G [0,T[ the support of f(t) is compact in M.^ , 
then this solution is unique. 

The proof of the theorem will be postponed to Section [5l It uses techniques 
of optimal transportation theory and extends the proof by Loeper in [^ for the 
uniqueness of solutions to the Vlasov-Poisson system. These techniques have been 
useful in the recent paper by Carrillo and Rosado [1] where uniqueness of solu- 
tions to several equations containing aggregation terms and aggregation/diffusion 
competition (e.g. swarming models, chemotaxis) was established. The main new 
difficulty in our proof is the vector potential A, which for a given / is defined by a 
non-linear integral equation. The existence of a solution to this integral equation. 
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the regularity of this solution, and the corresponding a-priori estimates used in our 
proof, require some elaborated work. 

The rest of the paper is organized as follows. In Section [5J we present some 
preliminary results on the scalar and vector potentials $ and A. We also discuss 
results on the linear Vlasov equation (|1.1|) and the associated characteristic system. 
In Section |3l we recall well-known results of optimal transportation that we use in 
Section|4]to establish some crucial estimates on the potentials, which are needed for 
the proof of Theorem [TJ Hereafter, L°°{X; Y) will denote the set of L°°- functions 
f : X ^ Y, and C^{X;Y) denotes the set of such functions of class C°° with 
compact support in X. By g G L°°{[0,T[,CbiX);Y), we mean that [0,T[3 t ^ 
g{t) e Cb{X;Y) belongs to i°°([0,T[). For a function M^ x rS 9 {x,p) ^ g{x,p), 
we denote by llyllioc (resp. ||g|JLoo) the L°°-norm of a; 1— > (7(2;,^) (resp. p ^^ g{x,p)), 
and by suppy the support of g. Also dxg (resp. dpg) denotes the matrix gradient 
of g with respect to x (resp. p). li D ~ [dij\ is a real matrix, we define \D\ — 

SUPij Mij|. 

2. Preliminaries 

Henceforth, we set the speed of light to c — 1. For simplicity and without loss 
of generality we shall omit the time dependence throughout this section, unless we 
specify otherwise. We first recall an estimate that will be used later on. For a proof 
see [I2I Lemma 2.7]. 

Lemma 1. For 1 < m < 3 let ro — 3/(3 — m) and r < ro < s. Set X = 

(1 — r/ro) / (1 — f/s)- Then there exists a constant C — C{m, r, s) > such that 
for any* £ L'^HL" {R^;R) 



ny)r^^ 



<C{m,r,s)\mlfmt^- 

In particular, C{m, 1, 00) = 3 (47r/rn)™' / (3 — m). 

The next lemma shows that the potentials $ and A of a weak solution to the 
RVD system are well-defined and satisfy the appropriate regularity. 

Lemma 2. Let f be given as in Definition]^ and define p according to ji.^[ ). // 
p € L°°{MP;M.) has compact support, then there exists bounded potentials $ and A 
that solve (Egj-fTT]]. Moreover, (a^*,^^.^) G Cfc(R^M x M^) ^nd therefore /V$ 
and /u^V^' are in ^-'^(R^; R'^), with va defined by il.l\ ). In addition, V • ^ = 
Trace(dj:A) — where 

dxA{x) = - / {w(g)WA - WA(8)w-H [3w®w- id] (da • w)}/(y,p)- 

^ JR3 Jr3 \y — X 

is the matrix gradient of A. 

Proof. The results corresponding to the scalar potential <I> are well-known, see for 
instance (T) Chapter 10]; therefore we center our attention on the vector potential 
A. First, we prove that the integral equation 

{2.1) A{t,x)^\ j [ [id + u;^u;]vAfit,y,p)-^^, va ^ ^^^ =, 



dpdy 
2 
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has a solution which satisfies the claimed regularity. The proof is similar to the 
one in [lH Lemma 7] for the classical solution of the RVD system. We reproduce 
it here for convenience. 

Let (7 > be a constant that may depend on /, to be fixed later on. Define the 

set 



-De 



A^Ck 



\A\\^^ < C} 



We show that there exists a fixed point A^o G Vq which solves (|2.ip . 

To start with, denote the kernel IC{x, y) = \y — x\~ [id + lu (E}uj] and let A G Vq. 
Note that \IC{x,y)\ < 2\y — x\~^ . Consider the mapping A n- r[^] defined by 
1 /■ /■ _ , p-A 



T[A]{x)^ 



IC{x,y)vAf{y,p)dpdy, va ^ 



l + \p~A\' 



We claim that T[A] G T>(j. Indeed, let {IC)ij(x,y) be the zj-entry of IC{x,y). For 
some Ml, U2 and ug, on the line segment between x and z, the mean value theorem 
implies 

(^).„- {x,y) - (/C).. (z,y) 
< 2 




Hence, since \va\ < 1, Lemma [T] implies 



\T[A]{x)-T[A]iz)\ < - 

< C\x 



\lC{x,y) -]C{z,y)\p{y)dy 
dy 



p{y): 



\y 



(2.2) < C{p)\x-z\. 

Thus, T[A\ is a continuous vector valued function. Also, it is a simple consequence 
of Lemma [1] that 

(2.3) \\T[A]\\^^ < 3(47r)i/3 ||p||^/^ \\p\\]!l ^ C. 

Therefore, T[A\ € Vq as claimed. 

Now, by virtue of the Schauder fixed point theorem [11' Theorem 3, Section 9.1], 
T has a fixed point A^o G T^q if T is a continuous mapping and the closure of the 
image of T is compact in 2?^. To show continuity, suppose that Ak ^f A in Pg. 

Since the mapping g ^^ v{g) = g [l + \g\ \ is C^, by Lemma [1] we have 



\T[Ak]{x)^T[A]{2 



9 
< C 



\vA, -VA\f{y,p) 
Ah^. 



dpdy 



\y-x\ 

< C{p)\\Ak 

To show that TVq C T>q is compact, we first notice that for A £ Vq, 

dy . „_ 1 



(2.4) 



\T[A]{x)\<\\p\\ 



supp 



p\x-y\ 



<C{p) 



1 
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Consider the sequence {Bn} C TV^ and let i? > be fixed. By ^^ and (|231) . 
the restriction 

{^f^)\{xeK^:\x\<R} 

is equicontinuous and bounded. Then, by Arzela-Ascoli's theorem and a standard 
diagonal argument we can find a subsequence {B^^} and a continuous, bounded 
limit vector field B such that {Bn^.} — > B uniformly on compact sets, and in 
particular pointwise. Clearly, ||i?||^cxj < C, and since {Bn^,} satisfies the estimate 
(|2.4p . so does B. We only need to show that the convergence {Bn^ } — )■ B is uniform. 
Indeed, let e > 0. Choose R > Q such that the right-hand side of (|2.4p is less than 
e/2 for |a;| > R. Then, for all k we have \Bnk{x) — B{x)\ < e for |a;| > R, and we 
can find a ko — ko{R, e) such that for all k > ko 

sup \Bn^{x) - B{x)\ < €. 
\x\<R 

This proves uniform convergence. Hence, all the hypotheses for the Schauder fixed 
point theorem are fulfilled, and thus T has a fixed point Aoo in Vq. 

Next, we show that Ac^ has the required regularity. To that end, define va^ 
and then ja^ according to (jl.ip and (II. 4p . respectively. Since \va^\ < 1, then 
\JA^\ < p. Thus, JA^ S L°°(R'^;M'') has compact support. On the other hand, 
1^(2;, y)\ < 2 \y — x\ , and since the imk-th entry of dxlC reads 

idxJC),,^,{x,y) ^ dx,{\y-x\-'[S.„. + ^W^]} 

= \y - x\^ [Stm^'' - (5fcmCj' - SikUj"^ + Sw'cj'^w™] , 

the kernel IC{x, y) satisfies the derivative estimate \dxK,{x, y)| < 6 |j/ — xp . Hence, 
we can use the standard theory for Poisson's equation to find that A^o G C^ (R^^; R'^) 
as claimed; see, for instance, Lemma 4.1] or [71 Theorem 10.2 (iii)]. The remain- 
ing assertions in Lemma [2] are easy to check. In particular, since |wAoo| £ 1 and 
p £ L°°{M.^;M.), we have by Lemma [1] that A^o and dxAo^ are bounded. This 
completes the proof of the lemma. D 

The potentials satisfy the following estimates: 

Lemma 3. Let f be given as in Definition [7] and define p according to {l.j^ . 
Suppose that $ and A are as in Lemma\^ If P ^ L°°(M'^;M), then 

(2.5) \\A\\^^ < C \\p\\%; ||p||i§ , lia.^ll^^^ < C ||p||^(^ \\pf/l . 

These estimates also hold for the scalar potential <&. Moreover, there exists a posi- 
tive constant C that depends on llpll^;^! and \\p\\i^oa such that , 
(2.6) 
\A{x) - A{z)\ + \dxA{x) - dxA{z)\ + |9,$(x) - 9,$(z)| <-C\x- z\ In \x-z\. 

for any {x, z) G M^ x R^ with \x - z\ < 1/2. 

Proof. These are standard results for the scalar potential $ which were al- 
ready used in [8] to prove the uniqueness of solutions to the Vlasov-Poisson system. 
Therefore, we only work here with the vector potential A. Following the notation 
in the proof of Lemma [H we have that \K,{x,y)\ < 2\y — x\~ and \dxlC{x,y)\ < 
6\y — x\~ . Then, since \va\ < 1 and thus \ja\ < p, the estimates in (|2.5p readily 
follow by Lemma [H 
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To prove (|2.6p we rely on a similar result discussed in jH Lemma 8.1] for the 2D 
Euler equation. To begin with, let h = \x — z\ < 1/2 and Br{x) be a ball of radius r 
centered at x. A lengthy but elementary computation shows that the kernel IC{x, y) 
and its derivative satisfy 



\IC{x,y)-ICiz,y)\<C 



y-x 



y- z 



\y-x\^ \y-zf 



1 



\y-x\ 



\y~z\ 



and 



\dylC{x,y) ~dylC{z,y)\ < C 



y-x 



y — z 



\y-x\^ \y-z\^ 



l2/-a;|^ \y-z\^ 



As a result, 



\Aix) - Aiz)\ + \d,A{x) - d^Aiz)\ 

< I {\IC{x,y)-IC{z,y)\ + \dylC{x,y)~dyIC{z,y)\)p{y)dy 



< 



y-x 



y - z 



\y-x\^ \y-z\^ 



y-x 



y - z 



\y-x\'^ \y-z\^ 



\y-x\^ \y-zf' 



= I + 11 + III + IV. 



\y-x\ \y-z\ 



p{y)dy 



In the remainder of the proof we shall only estimate I since it is slightly more 
involved than the other three integrals and they can all be estimated in the same 
fashion. To proceed, consider 



/ = 



+ 



lB2h{x) JB2{x)/B2h{x) 

= h+h + h- 
We estimate one integral at a time. 



V-B2(x) 



y-x 



y- z 



< 



IpIIl. 



dy 



B2h{x) \y-x\^ JB2h{x) \y- z\^ 



\y-xf' \y-zf 



dy 



p{y)dy 



j^ dr + J^ dr\ < C\\p\\ 



\x- z\. 



As for I2, let y G B2{x) / B2h{x) . The mean value theorem yields 



y — X 



y — z 



\y-x\^ \y-zf 



<C 



\x — z\ 

\y-u\^ 



for some u on the line segment between x and z. Then, since for some constant 
C > we have |y — a;| < C\y — u\, 

dy 



h < C||p||^^|x-z| / 
Je 

< C\\p\\^^\x-z 



B2(x)/B2h(x) \y-x\ 



dr 



2h T 



< -CIIpII^oo |a; - 2;|ln|a:: - z| . 
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To estimate /s, let y £ M.^ /B2{x). Then |y — xj > 2 and we use the mean value 
theorem and a standard estimate [9j Lemma 8.1] to find that 



y- 


-X y 






\y- 


|3 1 |3 

x\ \y-z\ 




< 


1 


1 




\y-x\ 


\y-x\ 


< 


\x-z\ 



\y-z 



y - 


- x 


1 


|2 


\y- 


-a,i 



^1 



\y-^ \y- 



\y-x\ \y- 




for some other u on the line segment between x and z. Hence, we have by Lemma 

m 



h < \x 



p{y)\ 



dy 



<C\\p\\ 



WPWh 



We gather these estimates and use the fact that |a; — z| < 1/2 to find that for some 
constant C[p) that depends on ||p||^i and ||p||^oo, 



(2.7) 



I < -C{p)\x-z\\n\x-z\. 



Thus, since the same rationale shows that (j2.7|) also holds for the integrals //, /// 
and IV, we conclude that 

\d^A{x) - d^A{z)\ < ^C{p) \x - z\ In \x - z\ , 

and the proof of the lemma is complete. D 

The next lemma characterizes the weak solutions of the RVD system via the 
associated characteristic system. We recall that the speed of light has been set to 
c= 1. 

Lemma 4. Let $ and A be given as in Lemma [^ Then, there exists a unique 
solution (X, P){s,t,x,p), < s <t <T, to the characteristic system 

(2.8) X = va{s,x,p), 

(2.9) p = [-V$ + i;^V^^] (s,2:,p) 

with {X,P)(t,t,x,p) — {x,p), associated to the (linear!) Vlasov equation il.l]) . 
Moreover, since the right-hand side of i2.8\) - (2^) is an incompressible vector field, 
the mapping {x,p) ^t- {X,P) is measure preserving. Conversely, the weak solution 
of il-l]) in the sense of U.5\) with /o G L"^ n L°°(M'';R), fa > 0, is uniquely 
determined by f{t,x,p) = fo{{X,P){Q,t,x,p)) on [0,T[/or all {x,y) £ B cR^, B 
a Borel set. 

1 /(J 

Proof. Set z = {x,p). Since the mapping g i— >■ v{g) = g \\ + Iffl") is C^, 

Lemma [3] guarantees that the vector field G = (ua, — V<i> + v\\/ A^^ on the right- 
hand side of the equations ()2.8|) - (|2.9|) is Lipschitz continuous in the momentum 
variable and Log-Lipschitz in space. This implies that there exists a unique solution 
[0,T[9 s ^ Z ^ {X,P){s, ■,-,■) to the characteristic system ([2^ - (f279l) and the 
characteristic flow Z{s,t,z) is Holder continuous with respect to z ^ Chapter 8]. 
Moreover, since by [Ml Lemma 1] we have in the classical sense 



(2.10) 



V,-G = V.^-VA + Vf 



-W<^ + v\WA' 



0, 
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the mapping M® 3 z i— > Z(-, -,2) -with inverse Z^^{s,t, •) — Z[t,s, ■)- is measure 
preserving [9l Chapter 8]. By (|2.10p . we can write the Vlasov equation in divergence 
form, i.e., 

(2.11) 5J(t,z) + V, •(G(i,z)/(i,z))=0. 

Then, for /o given as in the lemma, we have (as a corollary of Tj Theorem 8.2.1]) 
that the function /(i,z) = /o(Z(0,i, z)), t e [0,r[, z = {x,y) € B C K^ where B 
is a Borel set, with /I^^q = /o, is the unique solution to the equation (|2.1ip (and 
so dTTTj) ) in the sense of dHH). D 



Remark 1. Since by Lemma IH the weak solution / of (jl.ip in the sense of ()1.5|) 
satisfies / {t, Z{t, 0, z)) = /o(^) where z h- > Z(t, 0, z) is measure-preserving, then for 
all ip£ Co(IR^M), we have 

V{z)l{'t,z)dz^ I tp{Z{t,0,z))foiz)dz, 

which means that the map z n> Z(t,0, z) transports /o to /(i), or f{t,z) = 
Z{t,0,z)^fo as defined in Section |3] below; see p.2p . 

We conclude this section with an alternative formulation of the solution to the 
integral equation (|2.1I) . The following technical result will be useful. 



Lemma 5. Denote uj = (w^, w^,a;'^); Va; = (9i, 92,93) and let 5ik = 1 if i — k else 
5ik ^ 0. For any ip G C^(R3;K) anrf y £ R^^ w;e /lawe 

(2.12) / Aipix) [d,k - c^^^'] r^ = 2 / afc9,^(x)^^. 

Proof. Clearly, both integrals are well defined. Let the support of tp be contained 
in a ball centered at y with radius R > 0. Let < r < R. It is easy to check that 
9fcW* = —\y~x\~ [Sik - uj^u}''] and Aw* = — 2 |y - xp w* for |y - x| > r. To 
prove (|2.12p . we will show that 

(2.13) 2 / ^k^^(p(x)- — ^ = -lim / ALp(x)dkio'dx. 

Jr^ \y-x\ ^^Oj\y-x\>r 

Denote the integral on the right-hand side of (|2.13l) by /(r). In view of the compact 
support of tp, we can restrict the domain of integration of /(r) to r < \y — x\ < R. 
Then, integration by parts in dk and then twice for A yield 



/(r) = - / dk'p{x)Auj^dx + o{r) 

Jr<\y-x\<R 

= 2/ dk^{x)- ^+o(r), 

Jr<\y-x\<R \y~A 

where o(r) stands for all boundary terms at jx — y| ~ r. Note that the boundary 
terms at \y ~ x\— R vanish due to the compact support of Lp. It is not difficult to 
check that o(r) — ;► as r — >■ 0. Then, since we also have \y — xp w' = di\y — x\~ 
for \y — x\ > r, another integration by parts and a standard limiting process produce 
the identity (HHJ. D 
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Lemma 6. Let f be given as in Definitions^ Then, the solution A to the integral 
equation il.S]) for the vector potential satisfies 

AAix) = -47:jAix) - V (v • / m(2;)^^| 

I 7r3 \y-x\) 

in the sense of distributions, i.e., for any ip G C^(M'^;R), we have 
A{x)A(p{x)dx 



(2.14) 



= 47r / (p{x)JA{x)dx 



[ V^{x)\v- [ 

■/R3 L Jr- 



JAiv) 



dy 



\y-A 



dx. 



Proof. Since \va\ < 1, we have \ja\ < p and thus ja £ L°°(R'^;R'^) has compact 
support. Then, the integral in the curly brackets on the right-hand side of (j2.14p 
is of class C^(R"^;M'^). Therefore, the last term in (|2.14l) is well defined and so are 
the other terms of this equation. Now, substitution of (|1.3p into the left-hand side 
of (pji)) yields. 



A{x)/S.tf{x)dx 



Aip{x) 



[±d + uj®uj]iA{y) 



dy 



\y~x\ 



Aifi{x) 



1 . 

id — - (id — uj ®ijo) 



JAiy)-. 



dx 



dy 

\y-x\ 



dx 



(2.15) 
where, 

(2.16) 



h+h 



h 



Aip{x) < / JA{y) 



dy 



\y-x\ 



dx ~ All I (p{x)JA(x)dx. 



In the last step we have used that the integral in the curly brackets in (I2.16|) is a 
solution to the Poisson equation Ait = —Attja in the sense of distributions. As for 
the integral I2, we work per components: 



r 1 



Aip{x)\ / [S^k - uj'uj''] j\{y) 

UR3 

i / j^(y)(/ A^{x)[S,,-uj^uj''] 



dy 



\y-x\ 

dx 

\y-x\ 



dx 
dy 



JAiy)\ / dkd.,ip{x) 



dx 



\y-x\ 



dy 



(2.17) 



^^v{x) dk / j^Aiy) 



dy 



\y-^x\ 



dx, 



where in the second and third steps we have used Fubini's theorem and Lemma [5l 
respectively, and Fubini's theorem and integration by parts in the last step. Note 
that the resulting boundary term vanishes in view of the compact support of ip. 
The relation (pil| readily follows from ([^35)1 - (PTf)) . D 
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3. Tools from Optimal Transportation 

Denote by P2(M^ x M^) the set of probability densities f{x,p) on M^ x R^ with 
finite second moment, /ga^RS (kP + |p|^) f{x,p)dxdp < oo. The L^-Wasserstein 
distance between two densities fi{x,p) and f2{x,p) in P2(K'^ x M.^) is defined by 

Wiihj2) 

= infj / / \ix,p)-{y,q)\'^dj{{x,p),{y,q)); 7er(/ii,^2)} 

(3.f) = infj/ f \T{x,p)-ix,p)\^hix,p)dxdp; T#/i = ^j 

where d/ii = /i(a;,p)da;dp, d/i2 = f2{x,p)dxdp, r(/ii,/i2) denotes the set of all 
probability measures on M^ x MP with marginals /ii and /i2, and 7#/i = /2 means 
that 

(3.2) / / ip{x,p)f2{x,p)dxdp^ / / ip{T{x,p))fi{x,p)dxdp 

JR3 JR3 JR3 JR3 

for all test functions (p S Co(M'^ x M^). In [3j, Brenier proved that the minimization 
problem (I3.ip . the so-called Monge-Kantorovich problem, has a unique solution T, 
which is characterized /^i-a.e. by the gradient of a convex function (/> : M'^ x M'^ — > R, 
i.e., T is uniquely determined /xi-a.e. by T = V(f> with {V(f))^fi = /2 for some 
convex function 0. Note that in p.ip . the minimizers 7 and T are related by 
7 = (idR6,T)#/i. Moreover, if e [1,2] and 

(3.3) fe^Te^fi, Te = {2 - ejid^e + {6 - 1)T ^ V ({2 - 9)^-^ + (9 - l)c 

denotes McCann's interpolation ^lOj, then the curve [1,2] 3 9 t-^ fg e P2(R'^ x 
R'^) is the unique length minimizing geodesic joining /i to /2 in the Wasserstein 
space (P2(R3 x R^), VF2), in the sense that W2{hJ2) = W2{fije) + W2ifej2)- 
Furthermore, the interpolant fg satisfies the continuity equation in a weak sense, 

(3.4) dgfgix,p)+V^,p-iug{x,p)fgix,p))=0 V0G[1,2], 

where ug G i^^ (R^ x R'^; R'^) is the velocity field associated with the trajectory fg, 
i.e., 

(3.5) Ug {Tg{x,p)) = — '^—^ — = V(/i(x,p) - {x,p). 

Indeed, p.4p can be formally seen as follows. For any test function (/? £ Cg (R'^ xR'^), 
using p.2p with fg — {Tg)^fi and then (|3.5I) . we have: 

— / / ipix,p)fgix,p)dxdp 

W .03 ./n3 



Vip{Tg{x,p))ug {Tg{x,p)) fi{x,p)dxdp 

(3.6) =11 '^'Pix,p)ug{x,p)fg{x,p)dxdp 

ip{x,p)Vx.p ■ {fe{x,p)ug{x,p)) dxdp, 
where we use an integration by parts on the right-hand side integral in p.6p . 
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Combining p.ip - (|3.5|) . we have that 

W^(hj2) ^ f f \V^x,p)-{x,p)fhix,p)dxdp 
\ue{Te{x,p)) \'^fi{x,p)dxdp 

2 - 



(3.7) = \u0{x,p)\ fg{x,p)dxdp. 

Formula p.7p is commonly known as the Benamou-Brenier |2] characterization of 
the L^-Wassertein distance, namely, 

= mini / / / f{0,x,p)\u{e,x,p)fdxdpde; /(6l) G P2(M^ x R^)), 

where the minimum is taken over all absolutely continuous curves / : [1,2] 3 
e h-^ f{0) e P2(M^ X M3) satisfying the constraints /(I) = /i, /(2) = /a and 
Se/ + yx,p ■ [uf) = . For a development on this topic, we refer to Jg. 

In the next lemma, we collect some well-known results in optimal transport 
theory that will be needed later in the paper. 

Lemma 7. Let /i, /2 £ L°°(R3 x M^; r) ^g two probability densities in P2(M^ x R^) 
with compact supports. For any 9 G [1,2], define the interpolant fg as in i3.3\) . 
Then 

(i). For all d [1,2], fg has a compact support in Mp x M^, and 

(3.8) IIMIl.-, <niax{||/i||ic.^,||/2|U^~J. 

That is, [1, 2]3 9 ^ fg belongs to L°° ([1, 2], L^ n L°°{R^ x R^); R) . 
(ii). Moreover, [1,2] 3 ^^ fg is differentiate at every point 9 £ [1,2], and 
its derivative dgfg, defined in the weak sense by i3.6\) . satisfies 



(3.9) ||aeM|H-i(K^xR3) <max{||/i||L==^,||/2||Ls=J'/'W^2(/i,/2), V0e [1,2]. 
That is, [1,2] 3 9^ dgfg belongs to L°° ([1, 2], iJ-i(R3 x R3);R). 



Proof. The proofs of p.Sp and p.9p are done in [8J. Here we only show that fg 
has a compact support in R'^ x R'^. Indeed, assume that the support of fi, supp/^, 
is contained in the ball Bj^. centered at the origin with radius Ri for i — 1,2. By 
the definition p.3p of fg, Tg is the optimal map in W2{fi, fg)- Then for /ii-a.e., 
Tg is invertible and dx,pTg{x,p) is diagonalizable with positive eigenvalues (see [1] 
Thm 6.2.4 & Prop 6.2.12]). Moreover, the following Monge- Ampere equation holds 
for /xi-a.e. {x,p) e M^ x R^, 

fi{x,p) = fg {Tg{x,p)) detdx, pTg{x,p). 

It follows that for /ii-a.e. (x,p) e M^ x M^, if /e (T6)(a;,p)) 7^ 0, then fi{x,p) 7^ so 
that 

{{x,p) : fg{x,p) ^ 0} C Tg ({(a;,p) : f^{x,p) ^ 0}) 



i.e. supp/e C Te(B^J. But Te(5fl^) d {9 - 2)Br, + {0 ~ 1)Br, C Br where 
R = Ri + R2. Therefore supp/g C Br, i.e. /e has a compact support. D 
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4. Final Estimates 

We first note that if fo e L^ D L°°{R^ x R3;R) is sucli tliat /o > 0, and if / is 
a weak solution of the RVD system in the sense of Definition [U then f{t) > and 
II/(OIIli — II/o||li for all t e [0,r[, so that f{t) can be viewed as a probability 
density on R"^ x R'^ up to normalizing the L^-norm of /o to 1. Moreover, under the 
assumption that the support of f{t) is compact in M.'^ x M.'^ for all t e [0,r[, then 
fit) e P2(R3 X M3). 

For simplicity and without loss of generality, we shall omit the time dependence 
in / throughout this section. The next lemma gives estimates on the relative scalar 
and vector potentials of two solutions to the RVD system. 

Lemma 8. Let f\,h G P2 ni°°(R^ x R^; R). Define pi and jai according to \1.4^ , 
and let ^^ and Ai, i — 1,2, satisfy respectively the equations 

(4.1) A$,;(a;) = -47rp,(a;), lim $,(a;) = 

\x\—¥OQ 

and 

(4.2) AA,(x) = -4^j^,(x) - V (v • / JA,(y)^^| , , hm \A{x)\ = 0, 

in the sense of Ij.^p - JIT^ (see Lemma\^. Assume that for some R > 0, suppfi U 
suppf2 C Bfi X Bji. Then there exists a constant C > which depends on R, 
ll/illroo and ||/2||roo , such that 

(4.3) ||5,$i-9,$2lL2+||c',v4i-a,A2||i. <CT4^2(/i,/2), 
and, for i G {1, 2}, 

(4.4) / P^{x)\A^{x) ~ A2{xf dx <CWi{f^,f2). 

Proof. The estimate on the scalar potential, that is, the first term on the left- 
hand side of (|4.3p . is essentially proved in [8] under the weaker assumption of the 
boundedness of the charge density. Here we only prove the estimates on the vector 
potential, which are ()4.4p and the second term on the left-hand side of (|4.3p . Indeed, 
for i G {1,2}, set 

^•* \/i+i9.4.r 

By the assumptions on fi, we have that j^. G £^ fl L°°(R^;R^) and its support is 
included in Bn. Then, ja^ G L'^{M.^; R'^), and by the generalized theory of Poisson's 
equation O chap. 8], the Newtonian potential 



jAAy)y^eW'^Hl 
\y-x\ 



dy I ^ TT-i, 



which implies that 

v(v. / ^^^{y)-^KeH 

In view of Lemma [51 

(4.5) AA,(x) = -47rJA,(x)-v|v- / jA,(2/)^^j, i = l,2. 
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in the sense of distributions, and since the right-hand side of ()4.5p belongs to 
i/^^(M'^; R"^), we deduce by the generaUzed theory of Poisson's equation that Ai € 
W^''^{M.^;M.^). Taking the difference of the two equations in (14. 5p . we obtain 

~A{Ai~A2){x)^4tt{ja,~ja,){x) + v\w- [ {]a, - jA,){y)-, '^^ 



and since Ai ~ A2 G W^''^{M.^), then integration by parts against Ai — A2 yields 



/ \d^Ai-d^A2f dx^Air f {Ai-A2)- Ua^-] A,) dx+ f V -{A^- A2)V ■ I dx 

JK3 Jr3 Jt^3 

where 

Hx) = / (jAi -JA2){y) 



\y-x\' 

Notice that the boundary terms vanish. Indeed, standard arguments show that 
both dxAi{x) and dxl{x) have a decay Oda^P ), and Ai{x) has a decay 0{\x\~ ). 
Then, the products Aidxl{x) and Aid^Ai^x) have a decay 0(|x|~ ), which suffice 
to make the boundary terms equal to zero. By LemmaJU S/ ■ Ai — Q. Thus the last 
integral on the right-hand side vanishes as well. 
On the other hand, we write 



JAi-JA2 = / h{vAi-VA2)dp+ I VA^ifl- f2)dp 

= / fiivAi-VA2)dp- / VA2defed9dp, 

where fe is the interpolant (13. 3p between /i and /2. Inserting this identity into the 
above expression, we have 

/ \dxAi - dxA2f dx - At: f f h{Ai - A2) ■ (va, - VA2)dpdx 

= -47r / / [ {Ai~A2)-VA2defedpdxd9. 

Ji Jm Jr3 



Therefore, 



/ \dxAi - dxA2f dx ~ in / fi{Ai - A2) ■ {va, ~ VA2)dpdx 
= ^^111 fe^e ■ d. [(Ai - A2) ■ VA2] dpdxde 



(4.6) +47r/ / / f guP ■ dp [{Ai-A2)-VA2] dpdxde, 



where we use the continuity equation p.4p with the velocity field uo denoted by 
ug — (ug , Mg). Here Ug and Ug are the x and p-components of ug in R'^, respectively. 
We now estimate each of the integral terms in (14. 6p . First of all, consider the vector- 
valued function 

R3 9 z h^ v(z) = , e R^. 

Clearly, v € C^(K'^;R'^) and VAi — vigAi) where gAi{x,p) — p — Ai{x) and {x,p) £ 
Br X Br. It is easy to check that the derivative Dv{z) of v at any point z G K'^ is 



14 



REINEL SOSPEDRA-ALFONSO AND MARTIAL AGUEH 



given by the 3-by-3 real symmetric matrix 

1 



Dv{z) 



v/i + FF 



id 



1 



Then Dv{z) is real orthogonally diagonalizable [6j Theorem 2.5.6]. Moreover, if 
z = (zi, 22,^3) e K^, then detL'i;(2:) = (1 + \z\'^)^^/^ > 0, the determinant of 
the 1-by-l leading principal submatrix of Dv{z) is (1 + z| + z|)(l + |zp)~^/^ > 
(1 + |zp)~"^/^ > 0, and that of its 2-by-2 leading principal submatrix is (1 + z1){\ + 
|zp)^^ > (1 + |zp)~^ > 0, where the leading principal submatrices of Dv{z) are 
the upper-left square submatrices oi Dv{z). Therefore the matrix Dv{z) is positive 
definite [B] Theorem 7.2.5], and we have 



(4.7) 



Dviz)C-^>X\e V(z,e)e 



l^ X M^ 



where A > can be chosen as the lowest eigenvalue of Dv{z). By the mean value 
theorem, we have at every point {x,p) G Br x S^, 

{Ai-A2)-{vA^-VA2) = [v{gA^)-v{gA^)]-{Ai~A2) 

= Dv {gA2 + S{gAi ~ gA^)) {gAi - gA^)- (^1 - ^2) 
= -Dv (gA, + 5{gA, - gA,)) (^1 - A2) ■ {Ai - A2) 

for some 5 G (0,1). But since by Lemma [3l |gAi(a;,p)| < i^T^ for all {x,p) G BrxBr 
and for some Kr > 0, then 1(5^2 + H9Ai — gA2)) {x,p)\ < Kr uniformly on {x,p) G 
Bn X Br. Then using (|4.7p and the above identity, we have 



(^1 - A2) ■ {VA, - VA^) < -Cr\Ai 



A2? 



uniformly on Br x Br and for some constant Cr > 0. We then deduce that the 
left-hand side of (14.61) is bounded below as 



f \d^Ai - d^A2f dx ~ At: f f hiA^ - A2) ■ {va, - VA,)dpdx 

JR^ J Bn J Br, 



(4.8) 



>Cr [\\d,A^-d.,A2\ 



p\'\A,-A2) 



Ll 



for some other constant Cr > 0. On the other hand, by inserting the identities 
9^ [{Ai - A2) ■ waJ = VA^-{d^Ai-d.^A2) + {Ai-A2)-d^vA^ and9p [{Ai ~ A2) ■ va.^] = 



(^1 - A2 
side of 



dpVA2 into (|4.6p . it is easy to see that the expression on the right-hand 
is dominated by 



h + h = 



47r / / f fe \ue\ \d,Ai - d^A2\ dxdpdO 

fg\ug\\Ai-A2\\d^A2\dxdpde 



Att 



+ 47r 



fe\ue\ \Ai - A2I dxdpd6'. 



Since ||pe||i~ < (4/3)7ri?3 max{||/il|Loo, ||/2||loo} = K (see LemmalU), by Cauchy- 
Schwarz' inequality and Eq. p.7p . the integral /i can be estimated as 



(4.9) 



h < Ck Wd^Ai - 9,^2 11^2 W2{fij2)- 
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Similarly, using Cauchy-Schwarz' inequality and Lemma [71 we have that 
h < 4^VF2(/i,/2) / dof f [ fe\Ai-A2\^dxdp) 

< ^TtK W2{flj2) ( f \Ai^A2\^dx] 

\JB2R / 

and we deduce by Poincare's inequality that 

(4.10) h < Ck 119.^41-9,^2 11^2 W^2(/i,/2). 

a; 

As for /2 7 Cauchy-Schwarz' incquahty yields 

h < ^7rW2{fij2)\\d,A2\\]!l f d^( f I M^i - ^2|'dxdp) . 

" Jl Vjr3 jr^ / 

Then we use the second estimate in (j2.5p and Poincare's inequality to get, as for 

(4.11) h < Ck ||9,Ai -a,^2|li2 W^2(/i,/2). 
Combining (|46|) - (|4.10p . we have 

(4.12) Wd^Ai ^ d,A2\\l. + p\'\Ai-A2)^ <C\\d^Ai-d,A2\\^.W2{hJ2). 

for some constant C > 0, which inipHes that 

(4.13) ||9.Ai -a,^2|li2 < Ciy2(/i,/2). 
Finally, inserting (|4.13p in the right hand side of (|4.12p . we obtain 

py\Ai-A2) <CW2{flj2). 

LI 
This completes the proof of the lemma. D 

5. Proof of Theorem [T] 

For a non-negative function /o £ i^ n L°°(M'^ x M'^;M'^) with compact support, 
let /i and /2 be two compactly supported weak solutions of the RVD system with 
the same Cauchy datum /q. Let ($i,Ai) be the potentials induced by /i, i = 
1,2, respectively. Denote z — {x,p) and let < s < i < T. To ease notation, 
write Zi{s, t) instead of Zi{s, t, z) for the solution of the characteristic system ()2.8p - 
(|2.9p associated to the Vlasov equation (|l.ip . Equivalently, Zi{s,t) with inverse 
Z^^{s,t) — Zi{t,s) is the characteristic flow associated to the solution fi{t) = 
Zi{t,s)^fi{s) of the Vlasov equation; see Lemma S] and Remark [1] In particular 
Zi{t,0) is the flow associated to fi{t) = Zi{t,0)^fo, that is, fi{t,z) — fo{Zi{0,t,z)) 
for all z e -B C M®, -B a Borel set. We further denote Zi{t) instead of Zi{t,0) and 
define the function 

(5.1) Q{t) = l f fo{z)\Zi{t)-Z2{t)fdz. 

We have W^ {fi{t) J2it)) < 2Q{t) because 7 = (Zi(i), Z2(i))_^ /o is admissible in 
((3T|) : here the function (Zi(t), Z2{t)) : K^ ^ K^ xM*^ is defined by {Zi{t), Z2{t)) (z) = 
{Zi{t, z), Z2{t, z)). Clearly, Q(0) = 0. Our goal is to show that Q{t) = for every 
t G [0,r[. If so, then W2{fi{t),f2{t)) = which implies that /i = /a on [0,r[xR6 
and therefore uniqueness. 
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Take the time derivative on both sides of (|5.ip . By Lemma |4] we have 



Q{t)^ I foiz) Z,{t)-Z2{t) ■ Z,{t)-Z2it) dz 

6 L J L 

/o(z) [Xi(t)-X2(t)] ■ [vA,{i.Zi{t))-VAAt,Z2{t))'\dz 

fo{z) [ Pi(t) - P2{t) ] • [ V$i(t,Xi(i)) - V$2(t,^2(i)) ] dz 

+ I fo{z)\Pi{t)-P2{t)] ■ \v\^VA\{t,Zi{t))^v\^VA'2{t,Z2{t))]dz 
=: hit) + hit) + hit). 

In [S], it is shown that for some constant C > that depends only on HpiH^^oo , 
(5.2) hit)<CQit)il-\nQ{t)), 

provided ||.^i(i) — Z2(t)||^oo < e^^. This is essentially the result in [S] which yields 
uniqueness of weak solutions of the Vlasov-Poisson system under the assumption 
that the charge density stays bounded. As for the RVD system, it remains to 
estimate h and h- To estimate h, recall that 

VAit, Xit), Pit)) - viPit) - Ait, Xit))). 

1/2 

where vig) = g {l + g'^) . Then, since g i-^ vig) is C^, we have 

\vAAt,Xiit),Piit))-VAAt,X2it),P2it))\ 

< Ci\Piit)-P2it)\ + \A,it,X,it))-A2it,X2it))\) 

< C( \Piit)^P2it)\ + \Xiit)-X2it)\ + \Aiit,Xiit))-A2it,Xiit))\ ), 

where C — C(||pi||^oo ). Note the use of the second estimate in (|2.5p for the last 
step. Thus, Cauchy-Schwarz' inequality yields 



(5.3) h < C [^Qit) + Q^/\t)T^/^t)] 

where, in view of p.2p and (|4.4p in Lemma [51 

Tit) = [ foiz)\A^it,X^it)) - A2it,X,it))f 

JR6 



dz 



(5.4) 

with 

h < CQit). 



Poix)\A,it,Xiit)) ~ A2it,Xiit))\^ dx 
piit,x)\Aiit,x) — A2it,x)\ dx 

< cwiih,h), 



with C = CiR, WMl^ ). Therefore, since Wiifiit)j2it)) < 2Q(i), we find that 
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The integral I3 can be estimated just as h and h- Indeed, we have 

/0(z)[Pl(t)-P2W 

v\^it,Zi{t))~v\^it,Z2it)) I VA{it,Xiit))dz 
fo{z)\Pi{t)-P2{t) 
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■v\^{t,Z2{t))[WA\(t,X,{t))-VAt,{t,X2{t)) 

In view of (|2.5p . the first integral on the right-hand side can be estimated exactly 
as /i . On the other hand, the second integral on the right-hand side is analogous to 
I2, with the vector potential instead of the scalar potential. Hence, since {v^^l < 1, 
we can use mutatis mutandi the arguments in |8^ and Lemmas [3] and [5] to estimate 
this integral as in (|5.2p . 

Then, we gather all previous estimates to find that for some constant C — 

c{R, wnh^ )>o 



(5.5) 



Q(i)<CQ(i)(l-lnQ(t)), 



whenever ||^i(t) — ^2(0IIl°° ^ c ^- This is a Gronwall's-type inequality which 
yields Q{t) = on [0,r[ and therefore uniqueness. D 
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